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Abstract 

Induced charge-density (ICD) oscillations at the Cu(lll) surface caused by an 
external impurity are studied within linear response theory. The calculation takes 
into account such properties of the Cu(lll) surface electronic structure as an en- 
ergy gap for three-dimensional (3D) bulk electrons and a s — pz surface state that 
forms two-dimensional (2D) electron system. It is demonstrated that the coexis- 
tence of these 2D and 3D electron systems has profound impact on the ICD in the 
surface region. In the case of a static impurity the characteristic ICD oscillations 
with the l/yO^ decay as a function of lateral distance, p, are established in both 
electron systems. For the impurity with a periodically time-varying potential, the 
novel dominant ICD oscillations which fall off like ~ 1/p are predicted. 

Key words: Surface electronic phenomena. Noble metals. Low index single metal 
surfaces, Electronic surface states, Adsorbates 



In the last years, scanning tunneling microscopy has become a powerful tool 
to study phenomena governed by screening at metal surfaces, i.e., charge 
rearrangement in response to the disturbance caused by impurities or de- 
fects [1,2,3,4,5,6,7,8,9,10]. Remarkable examples of that are the visualization 
of surface-state-originated Friedel oscillations of the induced charge-density 
(ICD) at metal surfaces [1,2] and the investigation of the indirect long-range 
interaction between adsorbates mediated by these oscillations [11]. It has been 
demonstrated [12,13] that at noble-metal (111) surfaces this interaction, whose 
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energy decays as with the lateral distance p between adsorbates, can lead 
to mutual redistribution of adsorbed atoms up to p ~ 10 nm. This interaction 
might also be responsible for the rearrangement of adsorbate species on the 
Si(lll)--\/3 X -\/3-Ag surface [14]. Moreover, very recently it has been reported 
on a realization of a self-assembled two-dimensional (2D) atomic structure due 
to this interaction [15]. 

In the theoretical studies of such phenomena, the scattering approach has 
found the wide application [12,16,17,18]. Within this approach, a s—pz surface 
state at the noble-metal (111) surfaces is considered to form a 2D free electron 
gas, ignoring the fine structure of the surface-state wave function and the fact 
that this 2D electron gas coexists with underlying three-dimensional (3D) 
electron system. It is well known that the 2D electron system responds to 
the introduction of an impurity by producing Friedel oscillations with the 
characteristic 1/p^ decay [19], whereas in the 3D electron gas the oscillations 
fall off as R being a distance from the impurity [20]. Nevertheless, up 

to date, the question of how the 2D and 3D systems respond simultaneously 
in the vicinity of a metal surface has not been addressed yet. An answer can 
be obtained, in principle, from ab initio calculations. But, at present, these 
evaluations are feasible only for systems with the lateral distances between 
adsorbates of the order of several A [21]. 

In this Letter we investigate the response of an electron system at a metal 
surface to the introduction of an impurity, considering the Cu(lll) surface as 
an example. We show that taking explicitly into account the realistic Cu(lll) 
surface band structure, which is characterized by the energy gap at the center 
of surface Brillouin zone (SBZ) and the s — Pz partly occupied surface state, 
is crucial for the description of the surface response. 

In order to evaluate the ICD, we adopt an approach based on linear response 
theory in which an external perturbation Vext(i"';t^) and the corresponding 
ICD nmd(r; u) are related by the equation 



Here x(r, r'; u) is the density-response function which is non-local and complex 
and contains information on electronic excitations. In time- dependent density 
functional theory x(r, r'; uj) satisfies the integral equation [22] 



with x°(r,r';ti;) being the response function of a noninteracting electron sys- 
tem, t>c(r — r') is the Coulomb potential, and -ft'a;c(r, r'; u) accounts for dynami- 
cal exchange-correlation effects. As we are mainly interested in the evaluation 
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Fig. 1. (a) Surface electronic structure of Cu(lll). The grey area represents the 
projected bulk electronic structure. Solid line shows the surface state dispersion 
with an effective mass of m* = 0.42me. (b) The phase space available for e — h 
excitations. In the grey area the surface state intraband e — h excitations are per- 
mitted. The intraband e — h bulk excitations are possible everywhere. The inter- 
band "surface state — bulk" transitions are forbidden in the area below the blue line. 
Lines A, B, and C are respectively described by equations co = vjPq^\ + g'^/2m*, 

io = vjPq\\ - g^/2m*, and u = -vjPq\\ + q'^^/2m* with vjP being the 2D (surface 
state) Fermi velocity. Red line shows the acoustic surface plasmon (ASP) dispersion 
[27] very close to the line A (solid red line corresponds to the well defined ASP, 
whereas dashed one indicates the region of its gradual degradation). 

of long-range charge density oscillations, we use the random phase approx- 
imation, i.e., K^c = 0. The inclusion of a non-zero K^^ mainly affects the 
amplitude of Friedel oscillations [23,24,25]. 

To describe the Cu(lll) surface, we employ a slab containing 81 atomic layers 
of Cu together with a vacuum region corresponding to 20 interlayer spacings. 
The substrate is described by a model one-dimensional potential of Ref. [26]. 
This potential reproduces the key ingredients of the Cu(lll) electronic struc- 
ture shown in Fig. 1(a), namely, the energy gap at the center of the SBZ and 
the s — Pz surface and image-potential states in it. Using translation invari- 
ance parallel to the surface, one can perform the 2D Fourier transform for all 
quantities in Eq. (1). The ICD has now the following form [23] 



oo 
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nind{q\\,z]u) = / dz'x{q\\,z,z']u)Vext{q\\,z]u), 
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where Vext{<l\\, z'; uj) is the 2D Fourier transform of an external potential, Jo is 
the Bessel function of order 0, and z (p) is a distance perpendicular (parallel) 
to the surface. x° is given by^ 



id 

(ik|| /i,k|j+q|| ~ /i,k| 



(27r)2 -Ei,k||+q|| - -Ej,k|| +uj + iri 



(5) 



where the sum runs over the bands i and j, /'s are the Fermi factors, -E'i,k||+q|| = 
Ei + (k|| + q||)^/2m*, -Ej.kn = £j + ^/2"^j- £i {(pi) represent the one-particle 
energies (wave functions) obtained with the model potential [26] . The experi- 
mental values of effective masses m* different from the free electron mass have 
been used in Eq. (5). The regions of possible electron-hole (e — h) excitations 
at Cu(lll) are shown in Fig. 1(b). In the calculation of x"? a finite value for 
the broadening parameter rj = 0.001 eV was introduced. 

We begin with the case of response to a static external potential of the form 

Kxt(r) = -Qe/|r-ro|, (6) 



caused by an impurity of charge Qe = 1 placed at tq = {r\\ = O.zq}.'^ 
Fig. 2(a) shows the obtained ICD multiplied by + (2 — zq)'^ as a function 
of z and p. One can see that in a horizontal plane the Friedel oscillations 
peculiar to 2D electron systems are nicely reproduced: the ICD amplitude 
shows the quadratic dependence on p in the lateral direction and the period of 
the oscillations is precisely determined by the Fermi wave vector of the surface 
state, kjP = 0.1165 a.u. In the direction perpendicular to the surface the ICD 
oscillations demonstrate the behavior more complicated than that expected 
for a 3D electron gas. Note also that any displacement of the impurity along 
the ^-direction in the vicinity of the crystal border changes only the amplitude 
of the oscillations but not their character. 

To reveal the importance of the 3D system for the ICD shown in Fig. 2(a), 
we consider a hypothetical case of a free standing surface state without bulk 
electrons (the case of a pure 2D electron gas). Fig. 2(b) shows the resulting 
ICD. As in Fig. 2(a), one can observe the 2D Friedel oscillations in the lateral 
direction, while in the perpendicular direction the ICD behavior is in clear 
correspondence with the location of surface-state charge-density presented in 
Fig. 2(c). However, there is a qualitative difference between the oscillations 
presented in Figs. 2(a) and 2(b): the ICD as a function of z in Fig. 2(a) 

^ Hartree atomic units are used throughout unless otherwise is stated. 
^ In all figures, the impurity was placed at zq = 1.97 a.u. 
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Fig. 2. (a) Map of nj^^x + — 2o)^) (in a.u.xlO^) induced by an impurity placed 
at the point shown by shaded circle as a function of z and p. Horizontal dashed line 
shows the crystal border. The solid (vacuum) is located for z < (2 > 0). The 
surface atomic layer is at z = —1.97 a.u. (b) The same as in (a) for a hypothetical 
case of a free standing surface state without bulk electrons. Additionally, in (a, b) 
Tiind X (p^ + (-2 — -zq)^) as a function of z is shown by white curves for the vertical cut 
p = 87 a.u. (c) The surface state charge density as a function of z. (d) Shows the 
same as in (a) for the jellium model. Note the absence of lateral ICD oscillations in 
this case. 

changes sign whereas in the pure 2D case it does not occur (see white curves 
in Figs. 2(a) and (b) which give the ICD along the vertical cut p = 87 a.u.). 
Similar behavior is observed for any p. Inspecting Fig. 2(a), one can see that 
the ~ decay takes place both in the 2D electron system and in the region 
of the 3D system adjacent to the crystal border. In the latter system for a 
given z the ICD oscillates with the opposite sign to that in the 2D one. To 
clarify the origin of this ICD behavior we separate contributions to xo Eq. 
(5) into three parts: "surface state — surface state" {s — s), "bulk — bulk" {b—b), 
and "surface state — bulk" [s — b) transitions 



where a stands for s — s, 6 — 6, or s — 6, and the summation is performing 
over these three components. As a result, the ICD of Eq. (4) can be rewritten 
as the sum of the partial ICD concerned with the contributions of the 2D and 
3D systems as well as their mixing 
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where 



+ 
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/ 



dz'x°'{q\\,z, z'] u){Ve^^t{q\\,z'] uj) 



(9) 



Here Vc{q\\,z\z") 



2£g 911k' z"\ jg ^Yi^Q 2D Fourier transform of the bare 



Coulomb interaction, each x° is defined through Eq. (2) by the corresponding 



Thus, the ICD shown in Fig. 2(a) comprises three interrelated components 
nf^^ which represent the response of the corresponding "subsystems" to both 
the external perturbation (6) and the electrostatic field created by the two 
other nf„^. Note that if we neglect the Xs-b term in Eq. (7) and approximate 
the surface-state wave function by the (5-function, we obtain the results of 
Ref. [28]. The main advantage of Eq. (9) over Eq. (4) is the possibility to solve 
it iteratively and to analyze step-by-step how one component affects the rest 
and vice versa. Actually, by setting n\~^ = nl~^ = 0, we find n1~^ shown in 
Fig. 2(b). It can be considered as a first iteration. Then, by substituting the 
obtained nf~j* into n^~j' with n^~^ = 0, we find the response of the 3D system 
to the external field and the field created by the ICD of Fig. 2(b), and so on. It 
explains why the ~ decay of the ICD peculiar to the 2D system takes also 
place in the 3D electron system adjacent to the crystal border. Additionally 
the screening by the 3D system of the complex perturbation mentioned above 
is the origin of an alternate behavior of the ICD as a function of z at a given 
p (see the white line in Fig. 2(a)). 

To emphasize the role of the realistic electronic structure in the surface re- 
sponse, we plot in Fig. 2(d) the ICD obtained within a jellium model for 
Vs = 2.67 commonly used for description of Cu sp-valence electrons. In this 
case the distinct behavior of riind along the crystal surface without any os- 
cillations is clearly seen. Instead, this model gives strong ICD oscillations in 
the direction perpendicular to the surface which are suppressed in the re- 
alistic consideration (Fig. 2(a)) due to complicated mutual influence of one 
"subsystem" on the others. ^ Also, at small p, this influence leads to the ICD 
penetrating into the bulk only along a straight line at a finite angle away 
from the surface normal. Similar effect has also been observed in the studies 
of long-lived adsorbate-induced states at metal surfaces [29,30]. 

Now we consider the case of the potential (6) varying in time with some small 
frequency Uq as cos(ci;o^)- As we deal with the linear system, the superposition 
principle would hold true. Therefore, once the response to the given potential 

^ More detailed analysis will be given elsewhere. 
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Fig. 3. nij(g||,2; = 0;wo) vs g|| for four values of wq- tir for ujq = 0.22 eV, 0.44 eV, 
and 1.0 eV are displaced vertically by 0.1, 0.2, and 0.3 a.u., respectively. The arrows 
I (j) indicate the singularities positions at qA{QB,c)- 

is known, the response to an arbitrary time-varying external potential can be 
evaluated. In contrast to the previous static case when the integration in Eq. 
(3) over g|| is performed along the line u = of Fig. 1(b), here the integration 
should be performed along the line u = uq. Moreover, in this case the ICD 
depends on time and has the following form 



riindip, z; t) =nR{p, z; uq) cos{uJot) + n/(p, z; uq) sin(u;o^) 
= A{p, z] uq) cos(cjot - 0{p, Z] Uq)) 

where riR^p, z\ cjq) and ^/(p, z\ c^q) are defined through Eqs. (3) and (4) by the 
real and imaginary parts of the response function x(Q'|h ^''^ '^o); respectively. 
The ICD amplitude is A = \Jn\ + ri] and the phase shift 9 is defined by 
tan(6') = ni/riR. Thus, a non- vanishing nj leads to a non-zero phase shift 
between the ICD and external perturbation oscillations (it means an energy 
absorption in the electron system). 

Fig. 1(b) shows, as an example, three lines corresponding to t^o = 0.22 eV, 
0.44 eV, and 1.0 eV. Contrary to the "static" ujq = Q line which has only one 
singularity point at g|| = 2/c|P ensuring the Friedel oscillations in riind with 
the 1/p^ decay, the cuq 7^ lines have up to four singularity points (labeled 
by qA,B,c) due to the singularities at lines A, B, and C [31], what makes the 
dynamic response more complex. Commonly, these singularities are similar 
to that at q\\ = 2k'jP in the static case. But, at the Cu(lll) surface the 
qa singularity corresponds to the location of the acoustic surface plasmon 
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Fig. 4. Map of (a) and n/ (b) induced by the time-varying external perturba- 
tion (6) with frequency ujq = 0.44 eV as a function of z and p and multiphed by 
+ {z — zq)'^. Note the enhanced scale comparing with the one of Fig. 2. Other 
notions are the same as in Fig. 2(a). 
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Fig. 5. nR{p,z = 0;wo) multiplied by p"^ vs p for ujq = 0.0 eV, 0.22 eV, 0.44 eV, 
and 1.0 eV. On the left of the vertical line p = 0, the corresponding curves in the 
jellium model do not show any oscillations. 

(ASP) (see Fig. 1(b)), whose origin and nature has recently been studied 
[27,28]. The ASP singularity, appearing due to the coexistence of 2D and 
3D electron systems [28], is much stronger than the other ones and leads 
to the ~ cos{qAp)/p -like oscillations, i.e., with the spatial decay inversely 



proportional to p only. To visualize the distinct nature of singularities at 
and qB,c^ the dependence of riR on g|| for 2 = is shown in Fig. 3. 

Fig. 4 represents the ICD components nn (ICD at the moment t = 2n7r/uJo, n 
is integer) and rij {t = {2n + l/2)7r/u;o) for ujq = 0.44 eV. One can clearly see 
that in this case the ICD along the surface is dominated by oscillations with 
qa wave vector corresponding to the ASP. Nevertheless, the ICD also contains 
contributions originated from the other two singularity points, qs and qc- In 
general, the structure of the resulting oscillations is quite complicated, but 
fitted rather well by three cos-like functions. To compare the decay law of 
oscillations, we show riR^p, z = 0; uq) for uq = 0.0 eV, 0.22 eV, 0.44 eV, and 1.0 
eV in Fig. 5. For small Uq, the oscillations in 2D with = uo/vjP have the 
~ cos{qAp)/p asymptotic behavior. The same kind of ICD oscillations arises 
in the 3D system in the proximity of the 2D electron one. As ujq increases, the 
ASP looses its strength due to interband transitions (see Figs. 1(b) and 3), 
and oscillations with gradually change their decay behavior from to p~^. 
Thus, for ujq = 0.22 eV the decay obeys the p~^ law, while it is proportional 
to p~^-^ for uo = 0.44 eV, and p"^-^ for uq = 1.0 eV. 

In conclusion, we nicely reproduce within linear response theory the Friedel 
oscillations of the ICD at the Cu(lll) surface caused by a static impurity, 
taking explicitly into account the realistic surface band structure. It is shown 
that the coexistence of the 2D electron system with the 3D one has profound 
impact on the screening properties at the metal surface. Thus, additionally 
to the ~ cos{k'jP p) / p"^ ICD oscillations in the 2D system along the surface, 
the same kind of oscillations occurs in the 3D electron system in the region 
adjacent to the crystal border. In the case of a time- varying potential, dom- 
inant cos(u;o/'*^F^ ' p)/p "like oscillations at the surface in both the 2D and 
3D electron systems are predicted for low frequencies. We expect that these 
oscillations can lead to a longer-range indirect interaction between atoms and 
molecules at metal surfaces than it is thought nowadays. 
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